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ABSTRACT 

The production of massive quarks and leptons in e + e~ collisions is studied using ex- 
act helicity amplitudes. Total cross sections as a function of y cu t, in both the JADE 
and the kr algorithms, are presented and compared with massless results. Some 
invariant mass distributions are examined in relation to Higgs detection. Compact 
expressions for the helicity amplitudes are given. 
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1. Introduction 



Massive particles are abuntantly produced at e + e~ colliders. Often they are as- 
sociated with other, massless, particles in rather complicated final states. Whether 
masses are important, or whether they can be ignored, depends from the center 
of mass energy and also from the region in phase space which is under study. For 
instance, the electron mass can be ignored in most cases at LEP energies, but for its 
role in regulating collinear divergencies in reactions like Bhabha scattering at small 
angle. On the contrary, the top mass is so large that it must be included in all cases 
even at supercollider energies. In a recent paper [1] we have shown that at the Z° 
peak, when masses are properly taken into account, cross sections involving b or, 
to a lesser extent, c quarks differ significantly from the corresponding predictions 
obtained when masses are neglected. It is the purpose of this paper to extend these 
results in several directions. First of all we consider both LEP I and LEP II energies. 
At higher center of mass energy, mass effects are expected to be smaller, but only 
an explicit calculation can establish to which degree this is correct. On the other 
hand, one of the main task of LEP II will be to search for the Higgs boson in the 
mass range 45 GeV < run < 80 GeV [2]. Such a Higgs decays almost exclusively 
to bb and t + t~ and as a consequence it is important to determine as accurately as 
possible all reactions in which bb and t + t~ pairs are produced, which could provide 
a background to Higgs detection. Second, we extend the study of the production 
of up to four strongly interacting particles. Four-quark final states with quarks of 
different flavours are also examined. A number of jet-jet mass distributions are 
given, and the effects of different clustering algorithms are investigated. Third, we 
consider the production of two quarks and two leptons and of four leptons. This 
kind of events are relevant as a potential background to Higgs searches and provide 
a test of the Standard Model, though not at the level of precision which can be 
reached in more inclusive measurements. Fourth, we consider final states includ- 
ing photons which are actively studied experimentally in order to determine the 
electroweak couplings of the quarks, and as a mean to search for new phenomena. 
Finally we present the helicity amplitudes for all the reactions studied in this paper. 

We only consider processes without electrons in the final state, which proceed 
through e + e~ annihilation to photon or Z°. The two contribution are comparable 
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outside the Z° peak. Therefore both have been retained in all amplitudes. 

The powerful methods based on helicity amplitudes, which have enormously 
simplified the calculation of processes with several final state particles, can, in prin- 
ciple, handle massless and massive fermions with equal ease. In practice, however, 
the massless case has proven to be much simpler and a number of very elegant and 
physically relevant results have been obtained in this case. For a good introduc- 
tion and complete references, see [3] . The treatment of reactions involving massive 
fermions is still much less advanced. 

Several papers have studied three and four particle production at e + e~ col- 
liders [5-11] at the matrix element level but, to our knowledge, all of them lacked 
an essential ingredient, either considering only massless particles or neglecting the 
Z° contribution, for a consistent treatment of heavy particle production at LEP 
energies. 

The matrix element for all processes for which we give results has been com- 
puted at tree level following both the method of ref. [12,13] and that of ref.[14] as a 
check of the correctness of our results. In Appendix B we present the amplitudes in 
the formalism of ref. [12,13] which makes it simpler to give compact, easy to imple- 
ment formulae. In Appendix A, for the convenience of the reader, we briefly recall 
the method and the results of ref. [12,13] and we collect various formulae which are 
used in the analytic expressions. 

The amplitudes have been checked for gauge invariance, and for BRST invari- 
ance [15] in case of external Z°'s. In the appropriate limits our results reproduce 
those of ref. [4,9,11,26] after some misprints in the formulae of ref. [11] have been 
corrected. 

We have used M z = 91.1 GeV, T z = 2.5 GeV, sm 2 (6 w ) = .23, a s = .115, 
a em = 1/128, m T = 1.78 GeV, m c = 1.7 GeV and mj = 5. GeV in the numerical 
part of our work. 

In what follows we neglect all hadronization effects, and apply cuts at the 
partonic level. 

2. Jets 

All LEP experiments have performed a large number of QCD tests. For in- 
stance, mentioning only the measurements which are most likely to be sensitive to 
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masses, a s has been determined from jet rates [16] and flavour independence of the 
coupling has been verified [17]. Three [18] and four [19] jet distributions have been 
studied and compared with QCD predictions. The color factors, which determine 
the gauge group which is responsible for strong interactions, have been measured 
[20] . The possibility of tagging quark jets using the semileptonic decays of b and c 
quarks has been exploited, for example in studies of the differences between gluon 
and quark jets [21]. 

In the coming years, improvements in statistics, in secondary vertex recon- 
struction with silicon vertex detectors and in particle identification will allow much 
more detailed studies of heavy quark production at LEP. 

The experimental definition of a jet is based on a clustering procedure. The 
two most widely used schemes are the JADE algorithm [22] based on the variable 

v5 = 2^(i-cos<9 y ) (i) 

and the kx or Durham algorithm [23] which makes use of a new clustering variable 

rp min(R 2 R 2 ) 
Vl = 2 y^(l-cos^). (2) 

In fig.l we present the cross sections for e + e~ — > qqg and e + e~ — ► qqgg with q = d,b 
as a function of y cut for both definitions of y at LEP I. For small yj ut the cross section 
for bbg is almost 20% smaller than for ddg. As expected the ratio becomes closer to 
one for larger y cut , but for y cut as large as .2, still R^ d = a (bbg) /a (ddg) < .96 in both 
schemes. The cross sections for the same reactions at y/s = 200 GeV are given in 
fig. 2. The differences between the massive and the massless case are less important 
at higher energies, as expected, but still of the order of several percent. It has to be 
noticed that the Durham scheme tends to enhance these differences. This feature 
has been found in all cross sections and distributions we have studied. Therefore, 
depending on the kind of analysis which is performed on the data sample, different 
schemes can be used in order to exalt or suppress mass effects. A typical example 
of this behaviour is shown in fig. 3a and 3b where the invariant mass distributions 
in three jet events are presented. The different shape of of b and d events is quite 
noticeable for the Durham algorithm, not only for the invariant mass of the qg pair, 
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which is expected because of the different quark energy threshold, but also for the 
invariant mass of the qq pair. The exact shape of these distributions depends on 
Dent but the different sensitivity to mass effects in the two schemes remains. 

The qqgg final state is the dominant contribution to the four-jet cross section. 
The process with two hard gluon emissions is more sensitive to masses than the 
process with only one emission. For instance R^ d = a(bbgg)/a(ddgg) = .7 for 
yj ut = 1.5 x 10- 3 . 

The various contributions to e + e~ — > qqq'q' are shown in fig. 4 (q = q') and fig. 5 
{.1 7^ l') as a function of y cu t- Though very small compared with the qqgg rates these 
cross sections could be quite interesting if heavy quarks can be tagged with high 
efficiency. For 10 6 Z°, corresponding to an integrated luminosity of about 30 pb~ x , 
one expects approximately 100 events with 46 quarks and 400 events with 2b and 
2c quarks at y^ ut = 1. x 10~ 2 . The differences between massive and massless results 
are large, and even the charm mass has a significant effect. It is interesting to notice 
that the cross section for bbuu is larger than the cross section for bbdd, contrary 
to the naive expectations based on the fact that down-type quarks couple more 
strongly to the Z° than up-type quarks. This is due to the interference between 
graphs 1 and 2 with graphs 3 and 4 of fig. 12a. 

From fig.l, 2, 4 and 5 it can be seen that, while for three parton processes the 
cross sections for y^ ut = 1. x 10~ 2 are approximately matched by the cross sections 
for y^ ut = 1.5 x 10~ 3 , in the four parton case this happens for yj ut = 3. x 10~ 3 . 

The analytic expression of the matrix element for e + e~ — > qqg obtained in 
ref.[4], neglecting the Z°, has been exploited by some experimental groups [17] to 
estimate mass corrections to jet rate predictions which are used for measuring a s . 
The theoretical three jet rate to order 0(a 2 s ) for d quarks is multiplied by the ratio 
of the bbg to ddg tree level cross sections, R^ d after full detector simulation. We 
have compared the ratio obtained taking into account all contributions with the 
ratio obtained neglecting the Z°. The latter is consistently about 1% lower than 
the former at all y cu t- A further comment is in order. The theoretical three jet rate 
is a combination of three and four parton processes, the four parton contribution 
being of the order of 20% in the Durham scheme at small y cut . From our figures it 
is clear than the mass correction factors for three and four parton reactions differ 
considerably and therefore the simple method adopted in order to evaluate mass 
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corrections might slightly underestimate them. 

3. Backgrounds to Higgs searches 

The main production mechanism [2] for the Higgs particle at energies between 
170 and 200 GeV is e+e" -> HZ. If the Higgs mass is between 40 and 80 GeV the 
cross section for this process is of the order of 1 pb. As already mentioned the Higgs 
decays predominantly to bb. For an integrated luminosity of 500 pb -1 at y/s = 200 
GeV and ran = 60 GeV this yields about 30 events for e + e~ — > HZ — > bbl + l~ 
(/ = e 7/ u), 80 events for e + e~ — > HZ — > bbvv and 340 events for e + e~ — > HZ — > 
bbjj. Four-jet events are a potential source of background to this last type of 
processes. If 6's can be efficiently tagged, the main background contribution is 
given by e + e~ — * bbgg. For this reason we have computed the invariant mass 
distributions of all particle pairs in the bbgg final state. These distributions are 
shown in figs. 6, 7 and 8 at y/s = 91.1, 170 and 200 GeV respectively, in the two 
schemes. The intermediate point has been studied as a realistic energy for the first 
phase of LEP II. We have required yfj > 1. x 10 -2 or yf^ > 1.5 x 10 -3 for all pairs 
and that the angle of each particle with the beam satisfies | cos(#) |< .95. Again 
one notices the differences between the two schemes, particularly for small invariant 
masses. At LEP II energies the bb mass peaks well above 100 GeV while the gg mass 
clusters at small invariant masses. The cross section at y/s = 200 GeV, for 40 GeV 
< m-bb < 80 GeV is only about 5 x 10~ 2 pb which corresponds to approximately 25 
events. Moreover since the two gluons should fake an hadronic decay of the Z°, one 
could impose a cut on their invariant mass, decreasing this background drastically. 

The Higgs production cross section times the branching ratio to t + t~ at LEP 
II energies is about 80 fb. Hence it might be possible to detect the Higgs in this 
channel and to measure its coupling to the r. A comparison of this decay mode with 
the dominant one to bb would test the predicted proportionality of Higgs coupling 
to fermion masses. In fig. 9 we present the invariant mass distribution of the r + r~ 
pair in e + e~ — > r + r~qq at y/s = 170 GeV summed over five massless flavours. 
In order to select events which could fake an hadronic decay of the Z° we have 
required the invariant mass of the qq pair to be larger than 60 GeV. The cross 
section integrated for r + r~ masses between 40 and 70 GeV is less than 10 fb. 
Hence this background is rather small. On the other hand, with the mentioned cut, 
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a(e + e~ — > r + T~qq) = 86/6. This result is comparable to the Higgs production 
cross section in the t + t~ channel and it can be used as a reference point for the 
Higgs search. 

4. Jets plus Leptons and Four Leptons 

Four-fermion final states with at least one lepton pair are relatively easy to 
study experimentally, and in some cases they represent a potential background to 
Higgs production. Much interest has been spurred by the observation of an apparent 
excess in the t + t~X channel in ALEPH's 1989-1990 data [24]. This observation 
has not been confirmed by other collaborations [25]. 

Quite recently a new Montecarlo [26], which includes fermion masses and a 
complete treatment of 7 and Z° contributions as is appropriate outside the Z° 
peak, and combines them with initial and final state radiation and with a careful 
mapping of the many peaks in the matrix element due to collinear configurations, 
has been presented. We have computed the cross section for e + e~ — > / u + / u~r + r~ 
and e + e~ — > t + t~t + t~ , without initial state radiation, with the parameters used 
in ref.[26] and we obtain atot(e + e~ — > / u + / u~r + r~) = 640. ± 7. fb and atot(e + e~ — > 
t + t~t + t~) = 61.2 ± .3/6, in excellent agreement with that reference. Since we 
are primarily interested in massive particle production we have made no special 
effort to describe collinear configurations and have used the standard VEGAS [27] 
integration routine with the same parametrization over the whole phase space. Un- 
fortunately our straightforward approach does not allow us to obtain numerically 
stable results for e + e~ — > fx + \i~ fx + \i~ with the cuts of ref.[26]. Moreover we have 
not produced the additional diagrams needed to describe four-fermion production 
with final state electrons nor we have implemented initial state radiation, and this 
prevents us from a more thorough comparison between our results and those of 
ref.[26]. 

In table I we give the total cross sections for e + e~ — > t + t~t + t~ , e + e~ — > 
t + t~cc and e + e~ — > T + r~bb at ^/s = 91.1 and 200 GeV with the standard param- 
eters used in this paper, in particular with a em = 1/128 at all electroweak vertices. 
The expected number of events is very small, between 1 and 5 events per 10 6 Z° at 
LEP I. 



7 



5. Final states including a Photon or Z° 

Photons are the only particles which can be directly revealed and couple to 
the quarks in the early stages of the hadronization process. Final state radiation 
in multihadronic decays of the Z° can be used to test the electroweak couplings of 
up and down type quarks [28,29], combining the measurement of the radiation rate 
with the measured hadronic width of the Z°. 

In fig. 10 we give the cross sections for e + e~ — > qq^ and e + e~ — > qq^g, with 
q = (i, 6, at the Z° peak as a function of y cu t- The photon transverse momentum is 
required to be greater than 5 GeV. The difference between b and d quarks is again 
quite important, particularly for small y cu t- These differences are also visible in 
the photon energy spectrum given in fig. 11a and lib for the JADE and Durham 
scheme, respectively, even though the overall shape of the distributions are similar. 

In table II we present the total cross section for production of a Z° in association 
with a heavy fermion pair at y/s = 200 GeV. We implicitly assume that the Z° 
decays to a fermion pair different from the one that appears in the event. These 
processes are dominated by the diagrams in fig. 13b. Our results show that both 
intermediate photon and Z° give a sizable contribution. In fact if only diagrams 
with two Z°'s were important one would expect a(e + e~ — > Z(i + (x~) ps a(e + e~ — > 
Zt + t~); if, on the contrary, only photonic intermediate states were relevant one 
would get a(e + e~ — > Zee) ~ 3 x (2/3) 2 a(e + e~ — > Zt + t~). All processes produce 
rates of the order of a hundred events for 500 pb~ x of integrated luminosity. 

6. Conclusions 

We have computed the exact matrix elements, at tree level, for all processes at 
e + e~ colliders with three and four particles in the final state in which masses are 
relevant, with the only exception of reactions with final state electrons. We have 
studied these processes at LEP I and LEP II energies using both the JADE and the 
Durham clustering algorithm. 

Total cross sections involving b quarks are substantially smaller than the cor- 
responding ones for d quarks, particularly for small y cut . The effect increases with 
the number of jets and (obviously) with the number of massive particles. The dif- 
ferences between the cross sections or distributions for massive quarks and those 
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with massless ones are generally larger in the Durham algorithm. 

Four-jet events involving heavy particles and two-tau-two-jet events are a 
possible background to Higgs production. Therefore it is important to have exact 
expressions for all these reactions. We have studied several mass distributions in 
connection with the discovery of the Higgs and the measurement of its couplings. 
In the relevant mass regions these processes have rates much smaller than Higgs 
production. 

We have presented cross sections and energy distributions for one photon pro- 
duction in association with up to three jets, keeping quark masses into account. 
As in the case of purely hadronic final states, the overall shape of distributions is 
not strongly modified by mass effects, which however yield substantially lower rates 
than the corresponding reactions for massless particles. 

It will soon be possible, with improved statistics and tagging capabilities, to 
put to the test most of these predictions. 
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Appendix A: The helicity amplitudes method 



In this section, for completeness, we briefly recall the spinor techniques of 
ref. [12,13] which are used in our calculations. 

(i) Spinors. External fermions^ 1 ) of mass m and momentum p^ are described 
by spinors corresponding to states of definite helicity A, u(p, A) verifying the Dirac 
equation 



j/(u(p, A) = ±mu(p, A), 
u(p, X)^> = ±mu(p, A), 

and the spin sum relation 



(Al) 



J2 u(p,\)u(p,\) =jf±m, (A.2) 

where the sign +(— ) refers to a particle (antiparticle). 
One can choose two arbitrary vectors ko and k\ such that 

ko ■ ko = 0, k\ • ki = —1, ko ■ k\ = 0, (^--3) 
and express the spinors u(p, A) in terms of chiral ones w(ko, A) as 

u(p, A) = w(p, A) + /J,w(k , -A), (AA) 

where 

w(p, A) = jfw(k ,-X)/r], (A.5) 

and 



The spinors w(ko, A) satisfy 

(1) Unless stated otherwise, we shall use the term "fermion" and the symbol u 
for both particles and antiparticles. 
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w(k , X)w(k , A) = ^^/ o , (A.7) 

and therefore 

A=± ° 

The phase between chiral states is fixed by 

w(k ,\) = \jt w(h ,-\). (A.9) 

The freedom in choosing /c and hi provides a powerful tool for checking the 
correctness of any calculation. 

(ii) Polarization vectors for massless gauge bosons. External spin 1 massless 
gauge bosons of momentum p^ are described by polarization vectors corresponding 
to states of definite helicity A, A) satisfying 

e(p,X)-p = 0, e(p, A) • e(p, A) = 0, 



e"(p, -A) = e"*(p, A), e(p, A) • e(p, -A) = -1, (A10) 

and the spin sum relation (in the axial gauge) 

where g M is any four-vector not proportional to p M . 

Any objects £ M (p, A) obeying the relations (A.lO)-(A.ll) make an acceptable 
choice for the polarization vectors. For instance 

e»(p,\)=N[u(q,\)Yu(p,\)}, (A12) 
N being the normalization factor 



N=[4(q.p)]- 1 / 2 . 
11 



(A13) 



The existing freedom in choosing corresponds to fixing the gauge. The final 
results do not depend on the choice of q^. 

(iii) Polarization vectors for massive gauge bosons. For spin 1 massive gauge 
bosons there is an additional longitudinally polarized state satisfying 

e(p, 0) • p = 0, e(p, A) • e(p, 0) = 0, e(p, 0) • e(p, 0) = -1. (A.U) 
The spin sum becomes 

^(p, A)e"*(p, A) = -g^ + (A15) 

where m and p^ 1 are the gauge boson mass and momentum, respectively. 

For the polarization vectors of massive gauge bosons we cannot adopt the form 
(A. 12), since for a timelike momentum p^ 1 we cannot assign a definite helicity in a 
covariant way to the spinors u{p) and antispinors v(p). 

The simplest solution lies in noting that we are usually dealing with cross 
sections for unpolarized bosons, so we really have as the only requirement on the 
£ M (p, A) that their spin sum should be as in eq.(A.15). Any way by which we arrive 
at that expression gives us an acceptable choice for the polarization representation. 

Introducing the quantity 

a» = u(r 2 ,-)^u{r u -), (A16) 
where r%, r^ are two lightlike four-vectors satisfying 

r 2 = r 2 =0 , rf + r£=?A (A17) 

this can be used, after proper normalization, as the polarization vector belong- 
ing to p^ . In fact, if we replace the spin sum by an integral over the solid angle dQ 
of ri in the rest frame of the result is of the desired form: 

/ dn a^a v * = -nm 2 (-g^ + ^-). (A18) 
J 3 m A 
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This implies that we will obtain the correct result for the cross sections if we 
make the following replacements for on-shell massive bosons 



(iv) S and Z functions. Using the previous definitions one can compute 

S(X,Pi, P 2) = [u(pi, A)u(p2, -A)], (A.20) 

and 

Z(p 1 ,X 1 ;p 2 , A 2 ; p 3 , A 3 ; p 4 , A 4 ; c fl , c L ; c' R , c' L ) = 

[u(p u \i)T»u(p 2 , A 2 )] [u(p 3 , A 3 )r>(p 4 , A 4 )], (A21) 

where 

r (> = 7 M r (') 5 (^.22) 

and 

rn = ( c np fl + c n PL)j (A23 ) 

with 

ft = i^l, ft = (A24) 

the chiral projectors. 

The results (e 0123 = 1) are (see (A.6)) 

Q ,, n x Q (gi • fc o)(P2 ■ fci) - (pi • fci)(p 2 • fc ) + it^p^klp^pl 
b(+,Pl,P2)=^ , (A.Zb) 

S(- Pi,p 2 ) = S(+,P2,Pi)*, (A26) 

and 
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Z(p 1 ,+;p 2 , +;p 3 , +;P4, +; c R , c L ; c' R , c' L ) = 

-2[S(+,p 3 ,p 1 )S(-,p^,p2)c' R C R - ^i^mV^RCL - VlV2f-l>3lJ>4c' L C R ], 

Z(pi, +;p 2 , +;P3, +;P4, -; c R , c L ; c R , c' L ) = 

-2r7 2 Ci?[^( + ,p 4 ,Pl)^3CL - S(+,P3,Pl)^4<4], 

Z(pi, +; p 2 , +; P3, -; P4, +; cr, c l ; c^, c' l ) = 

-2r] 1 C R [S(-,p 2 ,P3)^4C L - S{-,P2,Pl)lA3C R \, 

Z(pi, +;p 2 , +;P3, -; P4, -; c fl , c L ; c^, c' L ) = 

-2[5 , ( + ,pi,p 4 )5'(-,P2,P3)c / L Ci ? - fJLi(l27l37l4</ L CL ~ VlV2fJ'3^c' R C R ], 

Z(pi, +;p 2 , -;P3, +;p4,+;c R ,c L ;c' R ,c' L ) = 

-2?74C^[5(+,p3,pi)^2Ci? - S(+,p 3 ,p 2 )^ic L ], (A.27) 

Z(pi, +;p 2 , -;P3, +;P4, -; c R , c L ; c' R , c' L ) = 0, 

Z(pi, +;p 2 , -;P3, -; P4, +; c fl , c L ; c^, c' L ) = 

£(pi, +;P2, -;P3, -;P4, -;c R ,c L ;c' R ,c' L ) = 

-2t]3c' l [S{+,P2,P4)h.iCl - S(+,pi : p±)ii 2 c R ]. 
The remaining Z functions can be obtained by exchanging + <-> — and R ^ L. 
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Appendix B: Matrix elements 



In this section we give the analytic formulae for the matrix elements. We define 
the propagator functions as 

1 1 

D iM = ^ D z(p) = sm 2 9wCOS 2 ew{p 2 - Mf + iM z T z ) ' {BA) 

and 



where M z and V z are, respectively, the mass and the width of the Z boson, 
Ow is the weak mixing angle and m/ is the fermion mass. 
We also define 

Ni = [4(qi'Pi)]- 1/2 , (B.3) 

where pi and ^ are the momentum and the auxiliary momentum of the z-th 
massless vector boson, respectively. 

Adopting for the polarization vectors of the gauge bosons the choices (A. 12) 
and (A. 16), and replacing any^( in the fermion propagator numerator with 

A=± 

one can express the Feynman amplitude T for a generic diagram as 

T = aCDM, (B.5) 

where a indicates the couplings, D the appropriate combination of boson 
and/or fermion propagators functions, C the eventual color matrix, and M a com- 
bination of Y and Z bilinear spinor functions. 

The Y functions are defined as 

Y(p 1 ,\ 1 ;p 2 ,\2;c R ,c L ) = [u(p 1 ,X 1 )(c R P R + c L P L )u(p 2 ,X 2 )]. (B.6) 
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Using (A. 4), (A. 5) and (A. 9) and computing the resulting traces one easily 

finds 



Y(pi,+;p2,+;c R ,c L ) = m x c R Vm 2 c L 



(B.7) 



Y{pi,+;p2,-;c R ,c L ) = c L S(+,p 1 ,p 2 ). 



(B.8) 



The remaining Y functions can be obtained by exchanging + <-> — and R <-> L. 

For all processes we will only report those spinor functions M which are not 
related by a trivial relabeling of momenta and helicities. We also adopt the symbol 
{A} to denote a set of helicities of all external particles in a given reaction and 
to indicate the usual sum over all possible helicity combinations. The expressions 
for the couplings cr and cl are given in table III. The amplitudes squared include 
the 1/nl factor for each n-uple of identical final state particles. Therefore the phase 
space integration must cover the whole space. 

I. Four fermion production: e~e + — > fff'f'. 



The electroweak production of two leptons and two quarks and of four leptons 



e-(pi, Ai) + e+(p 2 , A 2 ) -> r (p 3 , A 3 ) + £ + (p^ A 4 ) + /(p 5 , A 5 ) + /(p 6 , A 6 ), (B.9) 



is described in the case / = I (/ 7^ £) by the eight (first four) Feynman diagrams of 
fig. 12a plus the four (first two) of fig. 12b. The matrix element is given by 




e-e+ff 



(5.10) 



{A} J,m=l 




(S.ll) 



{A} J,m=l 
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with the amplitudes 

zT fc =e 4 J2 E D k( M k,vv ~ N k,vv) fe = 1, ...6(12), (S.12) 

where -Dfc denotes the propagator functions of the fc-th diagram. 

In formula (B.12) the spinor functions Nk,vv correspond to the p^p u term 
of the Z-propagator with momentum p, which is zero when acting on a massless 
fermion line. By definition Nk,vy = 0. For diagrams 5 through 8, 11 and 12 an 
additional minus sign has to be inserted because their are related to the remaining 
diagrams by the exchange of two identical fermions. 

We have 

M xyv ,= E z b5,A 5 ;P6,-A6;P3,A 3 ;Pi,A;4 v/ ,c{ v/ ;4^,ci v/ ) 

i=3,5,6 A=± 

x Z(pi , A; p 4 , -A 4 ; p 2 , -A 2 ; Pi , Ai ; c e Rv , c{ v ; ,c% v ), (B.13) 

Ni,w>= E E 

i=3,5,6 A=± 

[E E F b5,A 5 ;^,A / ;l,l)y(p :/ -,A / ;p 6 ,A 6 ,c^ / ,c{^)] 

j=5,6 A'=± 

X <E E [^(Pa.AaSPfc.V'a^jyCpfc.V'^^A;^,,^,)] 

fc=5,6 A"=± 

xZ(pi, X;p 4 , -X 4 ;P2, -A 2 ;pi, Ai; <4 v , c£ v ; c^,, cf, v )]}, (-B.14) 

Ai <-> -A 2 , A 3 <-> -A 4 , A 5 <-> -A 6 )]*, (S.15) 

Ai <-> -A 2 , A 3 <->• -A 4 , A 5 <-> -A 6 )]*. (-B.16) 
Then, the spinor functions corresponding to the diagrams of fig. 12b are 
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i=l,3,4 A=± 



x Z(p 5 , A 5 ; p 6 , -A 6 ; p 2 , -A 2 ; p», A; d R ^ , c£ v , ; c e Ry/ , cf, v , ) 

x Z(p 3 , A 3 ; p 4 , -A4; Pi, A; pi, Ai; , ci v ; , cf, v ), (5.17) 

where 

bi = -b 3 = -64 = -b 5 = -b 6 = -1. (5.18) 

Notice that 

Wfc+4(fc+8),vv' =0 fc = 1, 2(1, ...4), (5.19) 

since the Z-propagators are always connected to a massless electron line. 
Finally, in formulae (B.lO)-(B.ll) indicates the color factor 

C l = 1, or C q = 3. (5.20) 

The tree-level Feynman diagrams describing the reaction 
e~(pi, Ai) + e + (p 2 , A 2 ) -> q{p3, A3) + <z(p4, A 4 ) + g'(p 5 , A 5 ) + <f (p 6 , A 6 ) (5.21) 

are shown in fig. 12a (Only the first four contribute if q 7^ q'). The general form 
of the amplitude squared is 

WLw = I E E r W*Cb * (5-22) 

{A} Z,m=l 



l M lL" = 1^E E T^Ti^C^, (5.23) 

{A} l,m=l 

where C is a 8 x 8 matrix containing the color factors, and 



iT k =g 2 s e 2 J2 D kM k , v fe=l,...4(8), (5.24) 

v=~t,z 

are the individual amplitudes. 

The spinor functions and the propagators can be obtained from the correspond- 
ing expressions for the case e~e + — > ilff identifying V' — > g, £ — > g and / — > (/'. 
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Finally, the matrix of the color factors is 



C 



/ a 


a 


a 


a 


(3 


(3 


(3 




a 


a 


a 


a 


(3 


(3 


(3 


(3 


a 


a 


a 


a 


(3 


(3 


(3 


(3 


a 


a 


a 


a 


(3 


(3 


(3 


(3 






(3 


(3 


a 


a 


a 


a 


(3 




(3 


(3 


a 


a 


a 


a 


(3 




(3 


(3 


a 


a 


a 


a 


\(3 


(3 


(3 


(3 


a 


a 


a 


a / 



,ct = 2,P=- 



(5.25) 



II. Two fermion and one boson production: e e + — > ffV. 



For Z, 7 production in association with a // pair (/ ^ e) in 

e"(pi,Ai) +e + (p 2 ,\ 2 ) -> /(P3,A 3 ) + /(P4,A 4 ) + Z,j(p 5 , A 5 ), 



(5.26) 



the Feynman graphs are depicted in fig. 13a and 13b. For 7 production the 
amplitude squared is 



{A} Z,m=l 



(5.27) 



where 



-zT fc = e 3 ^ 5 fc M fcjV fc = 1, ...4. (5.28) 

V=7,Z 

Recalling that qi is the auxiliary quadrimomentum of the polarization vector 
corresponding to the massless boson of momentum pi, the y functions are 



M iy = N 5 ^ Z(p 5 ,X 5 ;q 5 , X 5 ;p 3 ,X 3 ;Pi,X; 1,1; c f R ^,c f L ^) 
i=3,5 A=± 

xZ(pi, X;p 4 , -X 4 ;P2, -A 2 ;pi, A x ; c^ v , c{ v ; c e Rv , c e Lv ), (5.29) 



M 2 ,^ = -[M ljV (pi ^p 2 ,P3 ^P4,P5 <-> 95; Ai <-> -A 2 , A 3 <-> -A 4 )]*, (5.30) 
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M 3 ,v = N 5 ^ b i XI Z (P3,A 3 ;p 4 , -A 4 ; Pi, A; pi, Xi, c f Rv , c f Lv ; c e Rv ,c e Lv ) 

i=2,5 A=± 

xZ(p 5 , A 5 ; q 5 , X 5 ;p 2 , -A 2 ;Pi, A; 1, 1; c R ,c e L ), (5.31) 



M 4) y = -[M 3) v(pi p 2 ,P3 <-> P4,P5 <-> <?5; Ai <-> -A 2 , A 3 <-> -A 4 ; &i <-> 6 2 )]*. 

(5.32) 

For Z production the amplitude squared has the form 

E^/«E«^ (-B.33) 

r \ l >0 ** l ~, 1 



{A} z " «,m=l 



where 



e 3 



-iTk = — ~ a t~ J2 D kM k ,v fc = l,...4. (5.34) 

V=j,Z 

The M/s 5 y functions can be obtained from eqs. (B.29) through (B.32) with 

the following substitutions: p$ —± r\, (75 — > r 2 , A5 — > — , iV~5 — > 1, — * and 

— > c Rz where r± and r 2 are the auxiliary momenta of the polarization vector 
for the Z. 

The color factor is 

C e = l, or C q = 3. (5.35) 



The Feynman graphs for 

e~(pi,Ai) + e+(p 2 , A 2 ) -> q(p 3 ,X 3 ) + q(p 4 , A 4 ) +^(p 5 , A 5 ) (5.36) 

are depicted in fig. 13a. The amplitude squared is 

{A} «,m=l 
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where 



-iT k =g s e 2 D kM k , v k = 1,2. (£.38) 

The spinor functions can be obtained from the corresponding expressions for 
e~e + — > qq'j substituting c q R ^,c q L ^ with c q R ,c q L . 
The color factor is 

C = 4. (£.39) 



III. Two fermion and two boson production: e e + — > ffVV'. 
The Feynman diagrams corresponding to 

e"(pi, Ai) + e + (p 2 , A 2 ) -»• g(p3, A 3 ) + g(p 4 , A 4 ) + y(p 5 , A 5 ) + g(p 6 , A 6 ) (£.40) 
are shown in fig. 14a and 14b. The amplitude squared for this process is 

WL = ^E £ (£.4i) 

{A} l,m=l 

with the amplitudes 

-iT k =g 2 s e 2 J2 D kM k , v fc = l,...8. (£.42) 

Recalling that ^ is the auxiliary quadrimomentum of the polarization vector 
corresponding to the massless boson of momentum p i: one gets for the spinor func- 
tions: 

M X y = N 5 N 6 Y £ £ £ ^(P5,A 5 ;g5,A5;P3,A 3 ;Pi,A;l,l;c^ fl ,4 g ) 

i=3,5 j=4,6 A=± A'=± 

xZ(p l ,X;p j ,X';p 2 , -X 2 ;Pi,Xi;c q Rv ,c q Lv ;c e Rv ,c e Lv ) 
xZ(p 6 , A 6 ; A 6 ;Pj, A';p 4 , -A 4 ; 1, 1; c Rg ,c q Lg ), (£.43) 
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M 3 ,v = N 5 N 6 Y Y Y z (P5,A5;Q5,A 5 ;P3,A 3 ;Pi,A;l,l;c^ g ,4 a ) 

i=3,5 i=3, 5, 6 A=± A'=± 



xZ(p 6 , A 6 ; q 6 , X 6 ;pi, X;pj, A'; 1, 1; c q Rg ,c q Lg ) 
xZ(pj,\';p4, -A 4 ;p 2 , -A 2 ; P i, Ai; c q Rv , c q Lv ; c e Rv , c e Lv ), 



(5.44) 



M 5 y = [M 3> v(pi <"> P2,P3 <"> P4,P5 <"> 96,95 <"> 

Ai <-> -A 2 , A 3 <-> -A 4 , A 5 <-> A 6 , )]*, 



(5.45) 



M 7y = N 5 N 6 Y Y Y { 2Z (Ps^ x s; Pi, X;p 5 ,X 5 ;q 5 ,X 5 ;c q Rg ,c q Lg ; 1,1) 

i=3,5,6 A=± A'=± 

xy(p 6 , A 6 ; p 5 , A'; 1, l)y(p 5 , A'; q 6 , A 6 ; 1, 1) 
+Z(p 5 , A 5 ; A 5 ;p 6 , A 6 ; q 6 , A 6 ; 1, 1; 1, 1) 
x [Y (p 3 , A 3 ; pe, V; 1, l)y(pe, A'; Pi , A; 1, 1) - y (p 3 , A 3 ; p 5 , A; 1, l)Y(p 5 , A; Pi , A; 1, 1)] 
-2Z(pe, Ae; ge, A 6 ;P3, A 3 ;p», A; 1, 1; c?j 9 , c| g ) 

xy(p 5 , A 5 ;p 6 , A'; 1, l)y(pe, A'; g 5 , A 5 ; 1, 1)} 
xZ(pi, A;p 4 , -A 4 ;p 2 , -A 2 ;pi, Ai; c q Rv ,c q Lv ;c e Rv ,c e Lv ), (BAG) 



M 8 ,v = -\M7,v{pi <-> P 2,P3 <-> P4; Ai <-> -A 2 , A 3 <-> -A 4 )]*. 
In this case, the matrix of the color factors is 



(5.47) 



C 



(1 


5 


7 


5 


7 


5 


e 




5 


7 




7 




7 


— e 


— e 


7 


8 


7 




7 


5 


e 


e 


5 


7 




7 




7 


— e 


— e 


7 


5 


7 


5 


7 




e 


e 


5 


7 




7 




7 


— e 


— e 


e 


— e 


e 


— e 


e 


— e 


c 


c 


\e 


— e 


e 


— e 


e 


— e 


c 


c/ 



,7=^,<5=-?,e = 6,C=12. 



(5.48) 
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The Feynman diagrams corresponding to the process 

e~(pi,Ai) + e + (p 2 ,A 2 ) -> <?(P3,A 3 ) + g(p4, A 4 ) + #(P5,A 5 ) + 7(p 6 ,A 6 ) (5.49) 
are depicted in fig. 14a and 14c. The amplitude squared is 

{A} «,m=l 

with the amplitudes 

-iT k = 9se 3 D k M ky fe=l,...10. (5.51) 

The spinor functions M\y through M§y can be obtained from the corre- 
sponding expressions for e~e + — > gggg substituting c|j , with , in the Z 
functions which have both ^6 and q& as arguments. 

The remaining spinor functions are 
M 7 ,v = N 5 N 6 ^ E ^(P5, A 5 ;g 5 , A 5 ;t>3, A 3 ;Pi,A; l,l;c|j a ,4 g ) 

1=3,5.7=2,6 A=±A'=± 

x Z(pi, A; p 4 , - A 4 ; , A'; pi , Ai ; c q Rv , c| v ; , c e Lv ) 
xZ(p 6 , A 6 ; q 6 , X 6 ;p 2 , -A 2 ;Pj, A'; 1, 1; cf^, 4 7 ), (5.52) 

M 8jV = -N 5 N 6 ^ 6 J Z (P5,A5;?5,A5;p 3 ,A 3 ;pi,A;l,l;c^ g ,c| g ) 

i=3,5j = l,6 A=±A'=± 

xZ(pj, A;p 4 , -A 4 ;p 2 , -A 2 ;pj, A'; , c£ v ; , c\ v ) 

xZ(p 6 , A 6 ; A 6 ;Pj, A';pi, Ai; 1, 1; c^, c^), (5.53) 

M 9jV = [M 8) v(pi ^P2,P3 ^p 4 ,Ai <-> -A 2 ,A 3 <-> -A 4 ;6i <-> 6 2 )]*, (5.54) 
Mio,v = [M 7)V (pi ^p 2 ,P3 ^P4, Ai <-> -A 2 ,A 3 <-> -A 4 ;6i <-> 6 2 )]*, (5.55) 
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where 



h=b 2 = -b 6 = -1. (£.56) 

Finally, the color factor is 

C = 4. (S.57) 
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TABLE CAPTIONS 

table I. Total cross sections for e + e~ — > t + t~t + t~ , e + e~ — > r + r~cc and e + e~ — > 
r + r"66 at = 91.1 and 200 GeV. Errors are as given by VEGAS [27]. 

table II. Total cross sections for e + e~ — > Z(i + (i~ , e + e~ — > Zr + r" , e + e~ — > Zee and 
e+e" -> at ^ = 200 GeV. Errors are as given by VEGAS [27]. 

table III. Right and left handed couplings of the fermions / = £, q to the gauge bosons 
1 ,Z,g. We have {ef,T(,gf) = (-1,-|,0) for / = e,fi,r; (e^T/,/) = 
(-i, -i, 1) for f = d,s,b and ( e /, t(, <//) = (|, 1, 1) for / = u , c. 
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Vtot 


(pb) 


channel 


LEP I (v/i = 91.1 GeV) 


LEP II (Vs = 200 GeV) 


T + T~T + T~ 


(876.4 ±4.3) x 10" 4 


(57.8 ±2.0) x 10" 4 


T + T~CC 


(2266.7 ±8.0) x 10" 4 


(269.5 ±3.5) x 10" 4 


T + r~bb 


(406.2 ± 1.7) x 10" 4 


(289.9 ±4.7) x 10" 4 



Table I 





<Jtot (pb) 


channel 


LEP II (Vs = 200 GeV) 




(341.1 ± 1.9) x 10" 3 


Z°t+t- 


(170.01 ±0.64) x lO- 3 


Z°cc 


(402.4 ± 1.0) x 10" 3 


Z°bb 


(378.19 ±0.71) x lO" 3 



Table II 
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7 


ry 
Zj 


9 


A 


ef 
ef 


—ef sin 2 6w 
T/ — ef sin 2 6*^/ 


9 f 
9 f 



Table III 
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FIGURE CAPTIONS 



fig.l. Cross sections for e + e~ — > bbg (continuous), e + e~ — > ddg (dashed), e + e~ — > 
bbgg (chain-dotted) and e + e~ — > <i<i<7<7 (dotted) as a function of y cu t for both 
definitions of y at = 91.1 GeV. 

fig. 2. Cross sections for e + e~ — > 66(7 (continuous), e + e~ — > cidg (dashed), e + e~ — > 
66(7(7 (chain-dotted) and e + e~ — > eieigg (dotted) as a function of y cu t for both 
definitions of y at ^/s = 200 GeV. 

fig. 3. Mass distributions for the qg pair and for the qq pair in e + e~ — > 66(7 (continu- 
ous) and e + e~ — > (iJ^ (dashed) events in the JADE (3a) and fcy (3b) schemes 
at y/s = 91.1 GeV. All particle pairs have yfj > 1. x 10" 2 (3a), yg > 1.5 x 10" 3 
(3b). 

fig. 4. Cross sections for e + e~ — > QgQg, q = d (dotted), -u (chain-dotted), c (dashed) 
and 6 (continuous) as a function of y cut for both definitions of y at y/s = 91.1 
GeV. 

fig. 5. Cross sections for e + e~ — > qqq'q' ', qq' = bd (dotted), du (chain-dotted), bu 
(dashed) and be (continuous) as a function of y cut for both definitions of y at 
v/i = 91.1 GeV. 

fig. 6. Two-jet invariant mass distributions in e + e~ — > 66(7(7 at y^s = 91.1 GeV. The 
pairs are 66 (continuous), bgi (dashed), 6(72 (chain-dotted) and gg (dotted) 
where gi ((72) is the most (least) energetic of the two gluons. All particle pairs 
have yf, > 1. x 10~ 2 (6a), yfj > 1.5 x 10~ 3 (6b). The angle of each particle 
with the beam satisfies | cos(0) |< .95. 

fig. 7. Invariant mass distributions as in fig. 6 at \fs = 170 GeV. 

fig. 8. Invariant mass distributions as in fig. 6 at \fs = 200 GeV. 
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fig.9. Invariant mass distributions of t + t~ pairs in e + e~ — > r + r~qq at -Js = 170 
GeV, summed over quark flavours. The invariant mass of the qq pair is required 
to be greater than 60 GeV in order to simulate a Z°. The quark masses in this 
case are neglected. 

fig. 10. Cross sections for e + e~ — > blr/ (continuous), e + e~ — > dd'y (dashed), e + e~ — > 
bb"fg (chain-dotted) and e + e~ — > dd^g (dotted) as a function of y cu t for both 
definitions of y at y^s = 91.1 GeV. All particle pairs have y > y cut . The photon 
transverse momentum is required to be greater than 5 GeV. 

fig. 11. Photon energy spectra in e + e~ — > blry (continuous), e + e~ — > dd"f (dashed), 
e + e~ — > bbyg (chain-dotted) and e + e~ — > dc^g (dotted) at y/s = 91.1 GeV. 
All particle pairs have y^ > 1. x 10" 2 (11a), yT. > 1.5 x 10" 3 (lib). The 
photon transverse momentum is required to be greater than 5 GeV. 

fig. 12. Feynman diagrams contributing in lowest order to e~e + — > qqq'q' (12a) and to 
e"e+ -> qqUoT e~e + -> M' (12a and 12b) with 7^ e. If the two fermion 
pairs are not equal only the first four diagrams in fig. 12a and the first two in 
fig. 12b contribute. A wavy line represents a photon or a Z° while a jagged 
line represents a gluon, a photon or a Z°. External lines are identified by their 
momentum as given in the text. 

fig. 13. Feynman diagrams contributing in lowest order to e~e + — > qqg (13a) and to 
e~e + — > (/(/F, F = 7, Z (13a and 13b). A wavy line represents a photon or a 
Z° while a jagged line represents a gluon, a photon or a Z°. External lines are 
identified by their momentum as given in the text. 

fig. 14. Feynman diagrams contributing in lowest order to e~e + — ■> <7<7<7<7 (14a and 14b) 
and to e~e + — > ^^(77 (14a and 14c). A wavy line represents a photon or a Z° 
while a jagged line represents a gluon or a photon. External lines are identified 
by their momentum as given in the text. 
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